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SYNOPSIS 


In the present work, solutions for the response of beams 
carrying concentrated mass, resting on soil medium have been 
presented idealizing the soil medium as Winkler model. Damping 
has also been included in these models. The beam with concentrated 
mass attached to it and resting on elastic foundation has been 
analysed for steady state dynamic force. The displacement 
magnitudes have been calculated for beams with different values of 
concentrated mass attached to it at different positions. It has 
been found, from the analysis, that as the concentrated mass 
increases the maximum displacement at the center of the beam 
decreases, for various positions of the concentrated mass. A 
general computer program has been presented for this beam. A study 
is also carried out for columns carrying concentrated mass at the 
centre of the column. A general computer program has also been 
given for these columns. 
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CHAPTER 1 


INTRODUCTION 

1 . 1 General 

The dynamic analysis of foundations finds many applications 
in modern technology such as power plants, machine tools, railway 
tracks, aircraft runways, buried pipes, rocket testing tracks etc. 
Such problems of foundation-soil interaction are generally 
analyzed by incorporating the reaction from the foundation into 
the response mechanism of the structure by idealizing the 
foundation by a suitable mathematical model. In majority of the 
cases, the response of the structure at the contact surface is of 
prime interest. Thus it would be of immense help in the analysis, 
if the foundation can be represented by a simple mathematical 
model with the Parameters of the model characterizing the true 
behavior of the system as closely as possible. To accomplish this 
objective, many foundation models, from the Winkler foundation to 
the elastic continuum idealization have been evolved. A 
comprehensive review of some of these foundation models is given 
by Kerr (1964). 

The dynamic loading effects on stationary structure have been 
the subject of great interest. A lot of work has been done on the 
behavior of beams subjected to steady state forces and transient 
loads. But for beams, resting on soil media the available 
literature is quite meagre. In technical applications it is not 
infrequent that a beam resting on soil media which is under 



vibration carries a lumped or concentrated mass somewhere along 
the length of the beam. A machine resting on a deep beam is an 
example. No solutions have been reported so far on such problems 
also, when the machine is eccentric to the beam. Thus it will be 
very useful in the analysis of the problems to consider the effect 
of the concentrated mass. 


1.2 Brief Review of Literature 

The earliest formulation of elastic foundation was due to 
Winkler who assumed that the foundation model consisted of closely 
spaced independent linear springs Kerr(1964). The Winkler model 
approximates the reaction of the continuum by a simple expression 
which results in a single ordinary differential equation governing 
the responses of the beam. Several other authors have evolved 
various Partly continuous models conceptually developed from the 
Winkler model, such as Pasternak, Filonenko Borodich, Weighardt, 
Hetenyi Kerr etc., [Kerr(1964) ] . The objectives in employing a 
model are mathematical simplicity and relatively accurate 
representation of the real system. Often, the visco-elastic 
model (Spring-Dashpot System) is used in dynamic analysis owing to 
mathematical simplicity, though the above mentioned models which 
are used for static loads, can also be used for dynamic analysis 
with appropriate inclusion of damping. In visco-elastic model, 
the soil medium is idealized in terms of physical constants- such 
as spring constant, damping coefficient and the equivalent mass of 
Soil participating in the vibration. The damping is attributed to 
radiation of energy due to the Propagation of the waves in the 




medium, and the energy dissipation through inter-granular 
friction. Through experimental evidence, Kausel and Rosset(1975) 
indicated that the energy dissipation takes place largely due to 
internal friction. It was also indicated that this energy is 
independent of frequency, about the spring constant, most of the 
methods presume that the spring constant is independent of motion. 
This assumption is not strictly valid [Whitman ( 1964) ] . However, so 
long as the motions are within certain limits, the behaviour of 
foundations may be taken as linear. 

Theoretical investigations on beam vibrations find its origin 
in Curvature-bending relationship derived by Bernoulli. 

For the transient loads, Stadler and Shreeves(1970) , Rades 
(1970) presented the analytical solutions for infinite and finite 
beam respectively for the dynamic response of the Euler-Bernoulli 
beam including the effect of damping, axial load or shear layer 
etc. , as given by equation 

+ T-2-X-+ Ky + c-?f + m-2-^ - q(x,t) (1.1) 
‘ 8x 2 3t St 2 

where T is axial load or shear layer coefficient. 

For beams an elastic foundation, some analytical results are 
available where the foundation moduli i are non-linear [Funstone 
and Hall (1967) ] . 

The problem of lateral vibration of a simply supported beam 
with concentrated mass attached at the centre was studied by 
Hoppman 11(1952), taking the Euler-Bernoulli beam. Chen(1963) 
presented a new formulation in which the dirac-delta function is 



used in differential equation. Frequency equations for normal 
modes of vibration were obtained by Baker(1964). Srinath and 
Das (1967) presented the analysis for simply supported beams having 
rotary inertia. Later Das and Deshmukh(1970) derived a general 
expression for Timoshenko beam with concentrated mass at the 
Centre. Jain(1977) derived a general expression for the response 
of beams on elastic foundations, with concentrated mass attached 
to it. 


1.3 Scope of the Present Investigation 


The present work is on the study of beams resting on soil 
media. The soil media is idealized as Winkler model with damping 
included in the model, for the analysis of Euler-Bernoulli beams. 
A beam with concentrated mass attached to it is analysed for 
transient vibrations, using finite element analysis. These 
solutions which can be used for any boundary conditions 
arechecked with the existing solutions for beams carrying a 
concentrated mass with end supports. Variation of magnitude of 
displacement along the length of the beam for various values of 
concentrated mass, mass position, is presented in this work. 
Comparisons have been made with some existing solutions. 

Also, the effect of concentrated mass on the axially loaded 
element is investigated. 

The general conclusions and scope for further work have been 


presented in chapter 4. 



Table 1.1 


Criteria for flexibility of beams 


S.No. Size of beam Criteria for Procedure 

distinction recommended 


1 . 

Long beam 

XL > 

0.5 

Winkler analysis 
assuming infinite beam 

2. 

Moderately long 

2.25 

< XL < 5.0 

Winkler analysis 

3. 

Short moderate 
beam 

0.8 • 

< XL < 2.25 

De beer method 

4. 

Short beam 

XL < 

0.8 

Treat as perfectly 


rigid. 
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CHAPTER - 2 


DYNAMIC RESPONSE OF BEAMS ON ELASTIC FOUNDATION CARRYING A 

CONCENTRATED MASS 


2.1 Introduction 

This chapter deals with the finite element solution of 
Euler-Bernoulli beam, carrying a concentrated mass attached to it, 
resting on elastic foundation. A machine resting on deep beam, a 
column resting on strip footing, a mechanical oscillator (to 
produce a high force oscillation) on footing, for the 
determination of soil properties are the few examples which 
simulate this problem. Also this chapter deals with, an axially 
loaded column on elastic foundation. 

In the literature, solutions for beams carrying concentrated 
mass and simply supported at the ends without elastic foundation 
are available [Hopmann(1952) , Chen(1963), Baker(1964), Srinath and 
Das (1967), Das and Deshmukh(1970) , Kameswara Rao and Das (1975)] 
and for beams carrying concentrated mass with elastic foundation 
has been given in [Jain(1977) ] . However the solution by 
Jain(1977), does not include the effects of eccentric concentrated 
mass. In this chapter a finite element formulation has been given 
for beams resting an elastic foundation, with Winkler model 
idealisation with damping included. Similarly, a finite element 
formulation for columns and carrying a concentrated mass has been 



given. 


2.2 General Theory 

The differential equation for transverse vibration of an 
Euler-Bernoulli beam on elastic foundation subjected to a dynamic 
load can be written as 


EI 

ax 4 


i ay . a 2 y 

^ - at ^ ^ 


= q (x, t) 


where 


E is Young’s Modulus of Elasticity, 
I is Moment of Inertia of the beam, 


k is Spring Constant of the soil medium. 


( 2 . 1 ) 


c is damping coefficient of the soil medium, 
pA £ is mass density of the soil medium, 

y is vertical displacement of the beam, 

EI is Modulus of rigidity of the beam, 
and q(x,t) is the dynamic forcing function. 

The above fourth order linear Partial differential equation 
is solved using the the finite element technique. As stated 
earlier, the soil is idealized by Winkler assumption with the 
effect of damping included. In the Winkler model it is assumed 
that the foundation applies only a reaction p(x) normal to the 
beam and that p(x) is directly proportional to the deflection 
i.e., p(x) = k y(x), where k is the Winkler foundation modulus in 



2.3 Finite element analysis 


For the finite element solution of Eq.2.1, the beam is 
divided into n number of finite elements as shown in Fig.2.1. A 
typical finite element is shown in Fig. 2. 2. 


Now, converting the Eq.2.1 into variational form, 
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(2.3) 


Similarly, 

l l 

J Ky Sy dx = 5 J y 2 dx 

O O 

1 2 1 2 
J pA ^ Sy dx = S J pA ^ y dx 
o c at o c at 


J C-|£-5y dx = sj C-|^-y dx 


(2.4) 


therefore the Eq.2.1 becomes 



i.e. , 

51 = 0 

where 



Now its solution is assumed as, 



= |N]{y> ne (2.7) 

where, ^ is the interpolating function of the element. 
Y is nodal parameter of the element 
r is the degree of freedom of the element 

also 

I I ne 

(y e ) = IN J {y> (2.8) 

(/ )" > iN"j W 


( 2 . 9 ) 


( 2 . 10 ) 


(y°) = [Nj <y> 

I ne 

(y- ) c = |Nj <y* > (2.11) 

The superscripts dashes , as shown above indicate the derivatives 
with respect to x, where as the superscripts dots represent the 
derivatives with respect to time. 


Now, the highest derivative of y in Eq.2.6 is of order two in 

the integral. Thus the interpolating function should have a 

(e) (e) 

compatibility of y and dy /dx i.e., the slope and deflection 
at the element interfaces should be continuous. The highest 
derivative in Eq.2.6 is of order three. Hence the interpolating 
function N^ should have a completeness of order three. 

y^ e ^ = a + bx + cx^ + dx^ (2.12) 

Upon applying the boundary conditions at element ends (y = y., y’ 

J 

= y* @ x = 0 and y = y, y* =y* @ x = 1), the Eq.2.12 becomes 



|N N N N I 

L 1 2 3 4 J 


f y 
y 


(2.13) 


where , 

N = 1-3 (x/1) 2 + 2 (x/1) 3 

N 2 = 1 ((x/1) - 2 (x/1) 2 + (x/1) 3 ) 

N 3 = 3 (x/1 ) 2 - 2 (x/1) 3 
and 

N 4 = 1((X/1) 3 - (X/1) 2 ) (2.14) 


Hence , 



I = s 1/2 El |yj (ne) <N"> L N ”j <y> (ne) dx 
o 

+ / 1/2 K l_yj (ne) {N>LNj{y> (ne) dx 

O 

* X PA C [Nj <y > (ne, lNj<y> l " e) <ix 
o 

+ X cl_Nj{y> (ne) LNJ <y> (ne) dx 
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- Lyj yj y k y k -1 


We have the condition. 



( 2 . 15 ) 
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(2.17) 
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Now this procedure is carried out for every element and the 
equation is to be assembled finally to bring it down to the form 

CM] {y} + [c] {y} + [K] {y} = {F} (2.18) 

There are two ways of solving Eq.2.18 : 

1. Direct Integration techniques and 

2 . Mode superposition technique 

Direct Integration Techniques : 

In this approach recurrence relations in time domain are used 
and depending on the type of relation various methods have been 
evolved as mentioned below. 

(a) Forward difference scheme, 

(b) Central difference scheme, 

(c) Backward difference scheme, 

(d) Houbolt method, 

(e) Wilson-8 method, and 

(f) Newmark' s method 

The methods (a) to (c) are explicit methods, and the methods (d) 
to (f) are implicit type. For these methods stability and 
convergence of the solution depends upon the time step chosen. 

In Mode super position, decoupled equations are obtained and 
then numerical integration is carried out. The Wilson-8 method is 


13 


used in the present work. Being an implicit method, this method 
is unconditionally stable, and hence, there is no critical time 
step and At can, in general, be selected many times larger than 
for central difference method. 


2.4 Wilson-© Method 

This method is essentially an extension of the linear 

acceleration method in which a linear variation of acceleration 

from time t to time t . is assumed. In the Wilson-© method the 
n n+l 

acceleration is assumed to be linear from time t to time t + ©At, 
where © £ 1, as shown in Fig. 2. 3. 


For © = 1.0 the method reduces to the linear acceleration 
scheme but in this case the method is only conditionally stable. 
For unconditional stability © a 1.37. For any time r so that 
0 s r s e At. We have from Fig. 2. 3 that 


y t+ r = y t + 


© At 


(y. 


t+© At 


y-> 


(2.19) 


Integrating Eq. (2.19) gives 


y t+ r - * t + y t T + 


2 © At 


(y 


t+0At 


- y-> 


( 2 . 20 ) 


and Integrating again gives 


y t+r = y t + ?t T + — y t T + {y 

C+T z z 2 6 © At 


t+©At 


YS 


( 2 . 21 ) 


or for the particular time r = 0 At, Eqs. 2.19 and 2.20 give 


t+0At 


= 5V + 


© At 


(y. 


t+©At + y t } 


( 2 . 22 ) 
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from which we can solve for y t+04t and J t+04t in terms of y t+0At = 


t+©At (e At) 


(y 


t+0At 




(0 At) 


V 2 y t 


(2.23) 


t+0At (6 At) 


(y 


t+0At 


- y t ) - 2 y t 


©At 


(2.24) 


To obtain the displacements, velocities and accelerations at 

time t+At, the equilibrium equations are considered at time t+0At. 

This requires projection of the applied load vector to time t + 0 

At, which is performed linearly as 

F = F + 0 (F . - F) (2.25) 

t+0At t t+At t 7 

and then Eq. (2.17) at t + 0At becomes 

I"> { y >t.8At + t C K?} t . eA t + [K1 M«® 4t ■ < F >t + 0At ,2 ' 26) 

Substitution of Eq.(2.23) and (2.24) into Eq.(2.26) gives the 

following system of simultaneous equations, which may be solved to 

s ive y t.eAt as : 


b [M] + b, [C] + [K] 

o 1 


M t . eAt - {F} t - « < F t , e4t - F t ) 

+ [M] (b o y t + b 2 y t + b 3 y t ) 

+ [C] (b t y t + b s y t + b s y t > (2.27) 

where 

b o = 6 /t 2 ; b t = 3/t ; b 2 = 2^ ; b 3 = 2 ; b 4 = 2 ; b g = x/2 
Substituting y t+e ^ t into Eq. (2.23) gives y t+ @^ t » which is then 
employed in Eqs.2.18 through 2.20, all evaluated at t = At. to 
give 
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y,.A, ■ Vy,.eA,- y,> * b ,«, + b ,y, 

S', .At ■ n - b , (y t . At ♦ y t ) 

s', .a, - y, + At y, +b , 0 ( y,.A, * 2 y,> 

where 


b = b / e ; b = -b / © ; b = 1 - 3/e 

6 o / 2 8 

b g = At/2 ; b iQ = At 2 / 10 


2.5 Effect of Concentrated Hass 

Now, for a beam carrying a concentrated mass and subjected to 
any generalized kind of transverse loading a finite element 
programme is developed. Fig. 2. 4 shows a beam carrying a; 
concentrated mass and resting on an elastic foundation, subjected 
to a steady state constant exciting force F o e . The inertia 
force is my. If this mass acts at element i, on nodes j, k, the 
elemental matrix due to this inertial force is £ m y{N} dx. For 
every element, the same procedure is carried out and is assembled 
finally. This assembled matrix is taken into mass matrix [m] on 
the left hand side, and the terms are added correspondingly. 

2.6 Finite Element Analysis of the Column 

From the Fig. 2. 5, the equilibrium of the element in the 
horizontal direction is, 


EA aU 

dx - m - — + q dx - fdx 

dX 2 

3t 2 

i.e, 


EA d * U 
dx. 2 

- m — ^ + q - f = 0 

3x 2 


0 


(2.30): 


A(x) is the area per unit length 
m(x) is mass per unit length 
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u(x,t) is displacement at any point x 

f is the coulomb friction due to soil per unit length 

Converting Eq.(2.30) into variational form 
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Now, its solution is assumed as, 


(2.31) 


(2.32) 


(2.33) 


(2.34) 


(2.35) 
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{u} e = L N N N Nl 

1 J u 1 2 3 r J 
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= LNJ{U}"' 


On completeness and compatiability considerations, u e : 

which upon applying boundary conditions at the element 

u at x = 0 and u = u at x = 1 ) , becomes 
J ^ 
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where , 


N = 1 - (x/1) and N z = x/1 


Hence , 
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Where g is concentrated load and £ is uni£orm iy distributed 
load. This de 3 ree coupled linear differential equation, can 
be solved on the same lines as has been discussed in the earlier 
section to obtain the dynamic response. 
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CHAPTER III 

RESULTS AND DISCUSSION 


3.1 Introduction 

This chapter deals with the numerical results obtained for 
beams and columns on elastic foundations carrying concentrated 
mass. In certain situations the length of the beam may be very 
long compared to its other dimensions viz., rocket test tracks. 
When the ends of the beam are far away from the load position, it 
may be assumed that the effect of load is negligible at the ends. 
This situation can be idealised as infinite beam. 

Response of beams and columns on elastic foundation carrying 
concentrated mass is considered in the present chapter. The 
solutions have been obtained by using Finite Element Method (FEM) . 
The load considered is a steady state constant exciting force. The 
programme developed can also take into account the earthquake kind 
of loads. Results obtained are compared with solutions by Jain [11] 
for steady state constant exciting force. 

The results are obtained for the concentrated mass in 
non-dimensional form at different positions along the length of 
the beam. The same has been done for different concentrated mass 
ratios (M/pA c L) . 

Results are obtained for columns without elastic foundations 
with concentrated mass at the centre of the column. 
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3.2 Finite Element Solution 


Consider a beam on an elastic foundation subjected to a 
steady state constant exciting force. It is required to determine 
the response of the beam. The partial differential equation of the 
beam is given by Eq.2.l. The corresponding equation for the finite 
element analysis is given by Eq. 2.17. Wilson-e technique is used 
for discretisation in time domain. 

A general purpose finite element programme is developed which 
can handle different boundary conditions, initial conditions, 
types of loadings etc. 

Also a programme is developed for columns which can handle 
different boundary conditions, initial conditions, types of 
loadings etc. 

The following data is used for beams and elastic foundations. 

6 2 

Young's Modulus, E = 2.1 X 10 kg/cm 

4 

Moment of Inertia of the beam, I = 1819.3 cm 
Spring constant of the 

2 

soil media, K = 525 kg/ cm 

Mass Density of the Soil media, pA c = 1.62 g/cm^ 

Damping Ratio, D =0.10 

Length of the beam, L = 1000 cm 

3.3 Verification of Finite Element Solution 

The finite element solution for an infinite beam resting on 
an elastic foundation and carrying a concentrated mass at the 
centre of the beam is verified with the already existing solution 
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by Jain [11], for a concentrated mass ratio of a = 0.34 and it can 
be seen from the Fig. 3.1 that the results are in good agreement 
with the results available. Also the results for beams without 
elastic foundation are verified as can be seen from Figs. 3.2, 3.3 
for simply supported and fixed beam cases respectively. 

The Response versus Time curve at the midpoint of the beam 
resting on elastic foundation, is plotted in Fig. 3. 4 for a 
concentrated mass ratio of a = 0.34, with a steady state exciting 
force of F q cos wt where F q = 210 kg. Similarly the Response versus 
Time curve for a concentrated mass ratio of a = 0.0 is plotted in 
Fig. 3. 5. It can be clear from these two figures that there is a 
decrease in the response with the presence of concentrated mass. 
Also Figs. 3. 6 and 3.7 show the variation of Response versus Time 
when the exciting forcing function is 274 cos wt, with and without 
concentrated masses respectively. It has been observed from these 
figures also that the response decreases in the presence of the 
concentrated mass . 

In Fig. 3. 8, the transient response is plotted for a beam on 
Elastic foundation with no damping in the system, for various 
points along the length of the beam. The dimensionless 
concentrated mass ratio a is 0.1 and the steady state constant 
exciting force is 274 cos cot at resonant frequency. The mass 
position is varied along the length of the beam, each curve 
representing a particular position from the left end of the beam 
as denoted in the Fig. 3. 8. It can be seen from this figure that 
when the mass position is 1/10 th of the beam from the left end 
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the response is higher than that of the no absorber case for 
points close to the centre of the beam and is almost negligible 
for points away from the midpoint of the beam. This is true for 
the mass position at 2/10 ths from the left end of the beam. But 
the response is almost same for mass positions x = 300 cms and 400 
cms from left end of the beam. For the mass position at the mid of 
the beam the response decreases and the rate of decrease is high 
as we move along points close to the centre of the beam. 

The same pattern has been observed for a = 0.2 (Fig. 3. 9) case 
also, but when the absorber is placed at x = 400 cms from the left 
end of the beam, there is a reduction in response, the. reduction 
being almost equal to the reduction observed in the mass position 
x = 500 cms. 

For a = 0.3 (Fig. 3. 10) it is observed that there is reduction 
in response even for mass positions from x = 200 cms onwards. The 
decrease in response increases as we move towards the centre of 
the beam. Similarly, the response curves are plotted for a = 0.4 
and a = 0.5 in Figs. 3. 11 and 3.12 respectively. 

Fig. 3. 13 shows the variation of transient response at various 
points along the length of the beam at various positions of the 
concentrated mass for oc = 0.10 with a steady state constant 
exciting force of 274 cos wt of resonant frequency. Here the 
nature of response reduction at various points is somewhat 
arbitrary, but when the concentrated mass is at the centre of the 
beam there is a reduction in the response and this reduction in 
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response increases as we move close to the centre of the beam. 
Similarly, for the concentrated mass ratios of a = 0.2, 0.3, 0.4, 
0.5, responses are plotted at various points along the length of 
the beam, in Figs. 3. 14 to 3.17, respectively. From these figures 
also it can be seen that the response is reduced when the mass 
position is at the centre of the beam and this response reduction 
increases as we move towards points close to to the centre of the 
beam . 

Fig. 3. 18 shows the variation of the transient response at the 
midpoint of the beam with different values of a, and with the 
absorber at the midpoint of the beam. It can be clear from this 
graph that there is a considerable amount of reduction in response 
as the value of a increases. For a. = 0.9, the vibration is almost 
reduced by 50% due to the presence of this mass at the centre of 
the beam. 

Response is also plotted for a column with one end free and 
one end fixed, without elastic foundation with concentrated mass 
ratios of a = 0.0 and a = 0.34 and it is seen from Fig. 3. 19 that 
the response is reduced at the free end of the beam. 
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CHAPTER IV 

CONCLUSIONS AND RECOMMENDATIONS FOR FURTHER RESEARCH 


4 . 1 Conclusions 

Based on the results obtained in this investigation the 
following conclusions can be drawn : 

1 . From the analytical results on the beam on elastic 

foundation carrying a concentrated mass on it and 
subjected to a steady state dynamic force, it can be 
seen that the effect of the concentrated mass is 
quite pronounced, for greater values of concentrated 
mass ratios. 

2 . The maximum displacement decreases very fast as the 
concentrated mass value increases on the beam. 

3 . The results obtained are in good agreement with the 
already existing results for beams carrying 
concentrated masses and resting on elastic foundation 
at the centre of the beam. 

4. When the mass position is at the centre of the beam, 
higher decrease in displacements are observed. 

5 . For columns without elastic foundation also the 
influence of the concentrated mass is observed on the 
displacement function of the column. 

6. The finite element method developed is very general 
and hence can be applied to a wide variety of beams 
on elastic foundation problems, just by changing 
input data in the programme. 
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4 . 2 Recommendations for Further Research 


1 . A spring mass and dashpot analysis can be done on 
both beams and columns on a elastic foundation case. 

2 . The transient response absorption is very important 
in earthquake kind of loading. Hence analysis for 
earthquake loads can be done with the absorber system 
used in the present work. 
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Fig. 2.1 Beam on elastic foundation divided into n 
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Fig. 2.2 Typical finite element of the beam on 
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Fig. 2.4 Beam on elastic foundation subjected to a steady state force 
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at resonant frequency L=1000.0cm 
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Response vs time at the center of the beam on an elastic 
foundation for a steady state constant exciting force of 210coswt 
with a concentrated mass ratio of a = 0.34. 
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Fig. 3.6 

Response vs time at the center of the beam on an elastic 
foundation for a steady state constant exciting force of 274coswt 
with concentrated mass ratio of a * 0.34. 
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Response vs time at the center of the beam on an elastic 
ramulu foundation for a steady state constant exciting force of 274coswt 

without concentrated mass. 
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